arXiv:1506.08303vl [nucl-th] 27 Jun 2015 


Null-test signal for T-invariance violation in pd scattering 

Yu.N. Uzikov“’^, A.A. Temerbayev'^ 

‘^Laboratory of Nuclear Problems, Joint Institute for Nuclear Research, Dubna, If 1980 Russia 
^ Department of Physics, Moscow State University, Moscow, 119991 Russia 
“L.N. Gumilyov Eurasian National University, Astana, 010008 Kazakhstan 

The integrated proton-deuteron cross section a for the case of the incident proton vector polariza¬ 
tion Py and tensor polarization P^z of the deuteron target provides a null test signal for time-reversal 
invariance violating but P-parity conserving (TVPC) effects. We study the null-test observable a 
within the Glauber theory of the double-polarized pd scattering. Full spin dependence of the ordi¬ 
nary strong pN scattering amplitudes and different types of the hypothetical TVPC pN-amplitudes 
are taken into account. We show that the contribution from the exchange of the lowest-mass meson 
allowed in the TVPC interaction, i.e. the p-meson, to the null-test observable a is zero. The axial hi 
meson exchange makes a non-zero contribution. We find that inclusion of the Coulomb interaction 
does not lead to divergence of the cross section a and calculate its energy dependence at the proton 
beam energy 100-1000 MeV. 
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I. INTRODUCTION 

CP violation (or T-reversal invariance violation under CPT symmetry) is required to explain the baryon asymmetry 
of the Universe [1]. In baryon systems violation of T-invariance has not been observed yet. CP violation established 
in physics of kaons and B-mesons leads to simultaneous CP and P-invariance violation. Under the assumption of 
CPT-invariance this implies existence of T-odd P-odd interactions. These effects are parametrized in the Standard 
Model by the CP violating phase of the Cabibbo-Kobayashi-Maskawa matrix. Another source for T-odd P-odd effects 
is the QCD 0-term, which can be related to electric dipole moments (EDM) of elementary particles and atoms in their 
ground states. 

On the contrary, time-reversal symmetry violating (T-odd) P-parity conserving (P-even) flavor conserving (TVPC) 
interactions do not arise on the fundamental level within the Standard Model, although they can be generated from 
the T-odd P-odd interaction by weak radiative P-parity non-conserving corrections. However in this case its intensity 
is too low [2, 3] to be observed in experiments at present. Thus, observation of the TVPC effects would be considered 
as indication of physics beyond the standard model. 

The existing experimental constraints on the TVPC effects in physics of nuclei are rather weak. So, the test 
of the detailed balance performed for the reactions Al{p,a)‘^‘^Mg and Al [4], and complemented by 

numerous statistical analyses of nuclear energy-level fluctuations leads to the ratio of T-odd to T-even matrix elements 
as ut < 2 X 10“^ [5]. Another type of experiment, i.e. polarized neutron transmission through a polarized ^®®Ho 
target gives ap < 7.1 x 10“"^ or pp < 5.9 x 10“^ [6]. Here Pp is the T-odd P-even coupling constant of the charged 
p-meson with the nucleon introduced in Ref. [7] to classify the TVPC interactions in terms of boson exchanges. 
Charge symmetry breaking determined as difference in scattering of polarized protons off unpolarized neutrons fm 
and polarized neutrons off unpolarized protons np gives ap < 8x 10“® (or pp < 6.7 x 10“^) [7]. One should add that 
indirect model-dependent estimation based on the existing constraints on EDM gives ap < 1.1 x 10“^ {pp < xl0“^) 
[8]. However, a more recent analysis showed [9] that EDM may arise via another scenario which suggests no significant 
constraints on the TVPC forces. 

The integrated cross section a will be measured at COSY [10] in double polarized pd scattering with a transverse 
polarized proton beam (pP) and a tensor polarized deuterium target [Pxz)- This observable provides a real null 
test of the TVPC forces [11]. This signal is not affected by the initial and hnal state interaction and therefore its 
observation would directly indicate time-invariance violation, as in case of the neutron EDM. The experiment [10] 
will be performed at a beam energy of 135 MeV. This energy choice was motivated by the theoretical analysis of the 
integrated pd cross section a performed in Ref. [12]. The aim of this experiment is to diminish the upper bound on 
the TVPC effects previously obtained in the n^®^Ho scattering [6] by one order of magnitude. 

The elastic channel and the deuteron breakup dp —> pnp were considered in Ref. [12] in the impulse approximation 
(single scattering mechanism) for estimation of a. In the present work we study the null-test observable d on the basis 
of the generalized optical theorem using the forward elastic pd scattering amplitude calculated within the Glauber 
theory. Both the single and double scattering mechanisms are considered. The spin-dependent Glauber formalism 
recently developed in Ref. [13] was applied in our previous work [14] to calculate spin observables of the elastic pd 
scattering using the strong (time invariance conserving and P-parity conserving) pN scattering amplitudes as input at 
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135 MeV. The obtained differential cross section, vector and tensor analyzing powers and spin-correlation parameters 
were found to be in reasonable agreement with the existing data [15, 16]. Here we generalize this formalism to allow for 
TVPC pN scattering amplitudes of several types. This generalized formalism is applied below to derive formulas for 
the null-test observable a and calculate its energy dependence. We show that within the single scattering mechanism 
this observable is zero in the Glauber theory (for any type of the TVPC pN interactions considered in the general 
case in Ref. [17]) and, consequently, focus on the double scattering mechanism. We investigate the contribution of 
several TVPC terms to the pN scattering amplitudes, in particular, the p-meson and axial hi- meson exchanges. In 
addition, we investigate the influence of the Coulomb interaction on the a cross section not considered in Ref. [12]. 

The paper is organized as follows. In Sect. II we consider the spin structure of the forward pd elastic scattering 
amplitude including the TVPC term and apply the generalized optical theorem to derive formulas for total spin- 
dependent cross sections in terms of the forward scattering invariant amplitudes. In Sect. Ill we construct the 
Glauber scattering operator taking into account full spin dependence of the elementary pTV-scattering amplitudes for 
strong and some types of TVPC interactions and S- and D- components of the deuteron wave function. Analytical 
expressions for the TVPC forward scattering amplitude g are derived for the double scattering mechanism with 
different TVPC terms. The influence of the Coulomb effects on the g amplitude is discussed in Sect.IV. Numerical 
results are shown in Sect.IV. 


II. FORWARD TRANSITION OPERATOR AND INTEGRATED CROSS SECTIONS 

Time-reversal symmetry conserving and P-parity conserving (TCPC or T-even P-even) interactions lead to the 
following transition amplitude of the elastic pd scattering at zero degree [18] 

e'p*ea = gi[ee'* - (me)(me'*)] -f 5 f 2 (me)(me'*) -f 
ig 3 {cr[e x e'*] — (crm)(m • [e x e'*])} -|- igi{crm){m ■ [e x e'*]), (1) 

where e (e') is the polarization vector of the initial (final) deuteron, m is the unit vector along the beam momentum, 
(T is the Pauli matrix, gi (i = 1,... ,4) are complex amplitudes. To the right-hand side of Eq.(l) one can add the 
TVPC (T-odd P-even) term in a very general form 

e'g*ta = 5{(cr • [m X e])(m • e'*) + (cr • [m x e'*])(m • e)}, (2) 

where g is the TVPC transition amplitude. To find the total spin dependent pd cross sections we use the generalized 
optical theorem [19] 




m- 


Tr {p^M{Q)) 
Trpi 


(3) 


where M(0) = 717(0)^*^^^ -I- M(O)^'^'^'^ is the transition operator from Eqs. (1) and (2) for the pd elastic scattering 
at zero angle 0 = 0, pi is the initial spin-density matrix, cr* is the total cross section corresponding to the density 
matrix pi. The transition operator is normalized according to the following relation with the differential cross section 
[13] 


— = -TrMM+. 
dt 6 

For the sum of Eqs. (1) and (2) one can write (see also [20]) 


(4) 


M{0)ap = giSap -k (32 - gi)mo,kp -k igscriCapi -k 7(^4 - 93 ) (JimimjCapj + g(Ti{f-zai'nipm^ -k e^pirriz^ma), (5) 

where (i = x, y, z) are the Pauli spin matrices, tap-y is the fully antisymmetric tensor, TOq (a = x, y, z) are the 
Cartesian components of the vector m. 

The initial state spin density matrix pi = pPp'^ is the product of the spin density matrices for the proton 

= ^(1+ P*’cr), (6) 

where is the polarization vector of the proton, and for the deuteron 

= 3 + 2^3P'j + gSjkPjk; 


(7) 
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here Sj is the spin-1 operator, pj and Pjk {j, k = x,y, z) are the vector and tensor polarizations of the deuteron, and 
Sjk = {SjSk + SkSj — ^Sjk) is the spin-tensor operator. Using Eqs. (3), (5), (6) and (7), one can find the total cross 
section of the pd scattering as 

<Jtot = al + alpP alipP ■ m)(p‘' • m) -h -f dpPP^^, (8) 

where (p'*) is the vector polarization of the initial proton (deuteron) and P^^ and Pxz are the tensor polarizations 
of the deuteron. The OZ axis is directed along the m, the OY axis is directed along the vector polarization of the 
proton beam pP and the OX axis is chosen to form the right-hand reference frame. The following equations were 
found in [21] for the TVPC terms 


o-Q = ^V^I'>n{2gi + g2), ct* = -A^/^^Img3, (9) 

0-2 = - 4 v^Jto (54 - 93), ^3 = 4 v^/m(gi - 32)- 

We find that the TVPC term g in the forward pd elastic scattering amplitude (2) leads to the following integrated 
cross section 


~ 2_ 

a = -A^/TTlm-g. (10) 

o 

In Eq. (8) the terms cr| with i = 0,1, 2,3 are non-zero only for the TCPC (T-even P-even) interactions and the last 
term a is non-zero if the TVPC (T-odd P-even) interactions occur. Thus, the term d constitutes the null-test signal 
for time-reversal invariance violating but P-parity conserving effects. This term can be measured in the transmission 
experiment [10] as a difference of counting rates for the cases with p^P^z > 0 and PyP^z < 0- 
We find the following matrix elements of the TVPC transition operator (2): 

< Ai' = i A' = = -1, A = 1 >= iV2g, (11) 

< /x' = i A' = = -i, A = 0 >= -i^g, (12) 

where p {p') and A (A') are spin projections of the initial (final) proton and deuteron on the beam direction, respec¬ 
tively. The diagonal matrix elements of the operator are zeros. For the operator the corresponding 

matrix elements in Eqs.(11) and (12) are identical and equal to V^g 3 - 


III. SPIN-DEPENDENT GLAUBER EORMALISM WITH STRONG AND TVPC INTERACTION 

A. Hadronic and Conlomb pN interaction 

Hadronic amplitudes of pN scattering are taken in a form of [13] 

Mn{p, q; cr, o-jv) = An + Cnctu + C^vcrjvh -|- BAr(crk)(cr^k) -|- (13) 

-f(GAr -f HN){eTq){(TNq) + {Gn - HN){eTh){eTNn), 

where q, k and n are defined as unit vectors along the vectors q = (p — p')) k = (p + p') and n = [k x q], respectively. 
Normalization of the amplitudes An, Bn, Cn, C'n, Gn, Hn is the same as in Ref. [13] 

J = \ttMnM+, (14) 

where da/dt is the differential cross section of the elastic pN scattering. The Glauber formalism for the pd elastic 
scattering accounting full spin dependence of the pN amplitudes (13) and S-and D-components of the deuteron wave 
function is given in Ref. [13]. The unpolarized differential cross section and analyzing powers of pd scattering calculated 
within this formalism are in reasonable agreement with existing data in the forward hemisphere at energies 250-1000 
MeV [13], [22]. Further development of this formalism was done in [14] to allow for calculation of spin correlation 
parameters and inclusion of the Coulomb interaction that is important at lower energies. 
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We include the Coulomb interaction by adding to the Glauber hadronic pd scattering amplitude the following pure 
Coulomb amplitude of the pd scattering 


MpM) = T^/w(q)'S'<i(q/2), 

txpp 

where is the antisymmetric Coulomb amplitude of the pp scattering [23]: 

4p(q) = fi^pp) - ^(1 + O- • 0-p)/(7!' - Opp) 


with 


f{dpp) — 


Avkpp sin^ Opp/2 


exp \—i— In sin + 2zxo 
V 2 


(15) 


(16) 


(17) 


here a is the fine structure constant, v {kpp) is the velocity (momentum) of the proton in the cms pp system, xo is 
the Coulomb phase. The momentum q transferred in the process pd —>■ pd is related to the pp scattering angle 9pp in 
the pp cms as q = 2kpp sin9pp/2. 

In Eq. (15) S'(i(q/2) is the elastic form factor of the deuteron which can be presented in the form [13] 

S'd(q/2) = So{q/2) - ^S' 2 ((?/ 2 )S'i 2 (q; (Tp, cr„). (18) 


Here 


5'i2(q; (Tp, cr„) = 3(crp • q)((T„ • q) - (Tp • (t„ (19) 

is the tensor operator, (T„((Tp) are the Pauli matrices acting on the spin states of the neutron and proton in the 
deuteron, q is the unit vector directed along the transferred momentum q. The form factors 5'o and S2 are related to 
the S- and D- components of the deuteron wave function [13]: 

So{q) = Sl°\q) + S^^\q), S 2 iq) = Si^\q) + Si^\q), ( 20 ) 

here 

poo poo 

So°\q)= dru'^{r)jo{qr),Sl^\q)= drw'^ {r)jo{qr), 

Jo Jo 

pOO i pOO 

S 2 \q) = ‘^J dru{r)w{r)j 2 iqr), S^^\q) =J drw'^ {r)j 2 {qr). ( 21 ) 


B. TVPC pN scattering amplitudes 

In general case TVPC NN interaction contains 18 different terms [17]. We consider here only the following terms 
of the t-matrix of the elastic pN scattering investigated in Ref. [12]: 

2 

tpN = fNia- ■ (TAr)(q • k)/mp + /iAr[((T • k){(TM ■ q) + {(Tn ■ k)((T • q) - ■ (T)(k • q)]/TOp + (22) 

+5Ar[(T X (Tat] • [q X k] /ml + g'j^ {cr - ctn) • i [q x k][T x TN]z/ml. 

Here cr (ctat) is the Pauli matrix acting on the spin state of the proton (nucleon N = p, n), x (x^r) is the isospin matrix 
acting on the isospin state of the proton (nucleon), q = p — pb In the framework of the phenomenological meson 
exchange interaction the term g'^ corresponds to the p—meson exchange, and the /lAr-term is caused by the axial 
meson exchange. As shown in Ref. [24], the contribution of the tt- and cr-meson exchanges to TVPC NN interactions 
is excluded. The TVPC NN interaction potential corresponding to /ii(1170), = 0~(1“' ) exchange in r-space 

has a form [25] 


Vh{r) 


GhGuml 

2X771^ 


Yl{x) [((Tlp)((T 2 f) 4- ((T 2 p)((Tlf)] 


(23) 
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where Gh {Gh) is the ordinary (TVPC) hNN coupling constant, nih is the /ii-meson mass, p = (pi— p2)/2, r = ri—r 2 , 
Yi{x) = (1 + x)e~^/, X = TTihr] cril2 and p^ are the spin-operator and momentum of z-th nucleon, respectively, and 
Vi is its r-coordinate (i = 1, 2). Making the Fourier-transformation of Eq. (23) we obtain the interaction potential in 
p-space and, therefore, find the factor hiq in Eq. (22) within the Born approximation for the tpN amplitude as the 
following 

2(^2 

hN = (24) 

where (j)h = Gh/Gh is the strength of the T-invariance violating potential of /ii-meson exchange relative to the 
T-conserving one, and PhNNio^) = (A^ — + q^) is the phenomenological monopole formfactor in the hNN 

vertex. Similarly proceeding from the TVPC p-meson exchange NN-potential in r-space [25-27] we find for the p'-term 
in Eq. (22) 


9n — 




2 w?p 


^pNN 


(q^), 


(25) 


where (ftp = gp/gp is the ratio of the TVPC pNN coupling constant Pp to the strong one Pp, nip is the mass of the 
p-meson, k is the anomalous magnetic moment of the nucleon, FpNN is the pNN vertex formfactor. 


C. The Glauber operator 


The Glauber operator of the elastic pd scattering in general case can be written as 

M(q,Q;S,<T) = jJJ e*Q‘-vI/+(r)OvI/,(r)d3r, 


(26) 


where 4'^ is the deuteron wave function, q is the transferred momentum, S = (cTp -|- crn)/2 is the spin-operator of the 
deuteron nucleons, (t/ 2 is the spin-operator of the incoming proton. We use the deuteron wave function generated by 
the NN interaction which conserves time reversal invariance and P-parity and has the following standard form 


4'd(r;(T„,crp) 


1 

■\/47rr 


u (r) -I- -^w (r) • 5*12 (f; cr„ 



(27) 


where lower index n (p) refers to the neutron (proton) of the deuteron target; u and w denote the S- and D- wave of 
the deuteron, respectively; the tensor operator S '12 (f;cr„,crp) is defined by Eq. (19). 

The operator for the single and double scattering mechanisms of pd scattering in general case (beyond the collinear 
kinematics) can be written as an expansion over the Pauli matrices (t„, cvp and in notations of Ref. [13] takes the 
form 


O (cr, (Tn, (Tp) =U {(t) + V„ (cr) ■ (Tn+Np (ct) ■ (Tp + Wij ( ct ) • {aniCTpj + UnjCTpi) , (28) 

here z, j = n,k are indices of the projections onto directions of three orthogonal vectors q, n, k introduced in Eq. 
(13). The operators C/, V, W act only on the spin state of the beam proton and do not depend on the spins and 
coordinates r of the target nucleons. When making the matrix element of the operator (28) over the deuteron wave 
functions (27) we obtain from Eq. (26) 


1 

7 i 


M(q, Q;S,cr) = 


g*Qr^+(r)o^.^(r)d3r 
1 


: USo + VS5'(°^ + {Si, Sj} - Wu] S^°'> - -^USuiQ; S, 8)^2 


Si2{Q;V,S)si^'> + V8WuSi2{Q;S,S)si^'> - ^5i2(Q;V,S),sf) - 

-lE.,5i2(Q;S,S),sf) 


( 2 ) 


-V2W,j [{5i, S,} Si2{Q; S, S) + Si2{Q; S, S) {S,, S^} 

+ J ^12 (f; <t„, <Tp) S^} Su ( 


5 ^ + 




( 29 ) 
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Here we use the notations V = Vp + V„ and = SiSj + SjSi; the form factors 

So{Q), S 2 {Q), Sq^\q), S 2 ^\q) and S^\q) are defined in Eqs. (20), (21); the tensor operators S'i 2 (Q;V,S) 

and S'i 2 (Q;S,s) are defined in Eq. (19). In Eq. (29) summation is performed over repeating indexes i,j. To make 
the integration over directions of the vector r in Eq. (29), we used the following relation [28] 

JJ dn^ exp {-iClr)Ti{r) = 47rj/(Qr)(-z)'T/(Q), (30) 

where ji{x) is the spherical Bessel function, r 2 (n) = (cr^ • n)(cr„ • h) — ^(cr^ • cr„), To(fi) = a-p ■ n, Q and f are 
unit vectors along n, Q and r, respectively. 

Eq. (29) is a generalization of Eq. (18) from Ref. [13]. The difference from Ref. [13] consists in two following 
respects. First, the operators U, V, Wij contain not only T-even P-even terms, but T-odd P-even terms as well. 
Second, we present in Eq. (29) all terms allowed within the Glauber theory, whereas in Ref. [13] small spin-dependent 
terms (of the order higher than two in definitions of Ref. [13]) were neglected. These terms are small at high energies 
about 1 GeV, however, may be important at lower energies ^ 100 MeV, corresponding the GOSY experiment [10]. 


D. Differential spin observables 

Vector analyzing powers Ay and spin correlation coefficients Cij, Cij^k of the elastic pd scattering are calculated as 

Al = TrMayM+/TrMM+, = TrMSyM+/TrMM+, 

C^z,y = TrMS^^ayM+/TrMM+, Cy^y = TrMSyayM+TrMM+, 

= TrMS^a^M+/TrMM+, = TrMS^a^M+ /TrMM+, (31) 

where M is the transition operator given by Eq. (26). Details of these calculations in terms of invariant amplitudes 
are described in Ref. [14]. 


IV. NULL-TEST SIGNAL OF TVPC INTERACTIONS 

Eq. (29) gives the single scattering pd amplitude if one put Q = q/2, where q is the momentum transferred in the 
pd scattering. Within the Glauber theory the amplitude of the single scattering mechanism is proportional to the 
on-shell tpAr(q) amplitude. At zero scattering angle the TVPC amplitude (22) vanishes, therefore, the corresponding 
pd-scattering amplitude and total cross section a of pd-scattering are equal to zero in the single scattering Glauber 
approximation. Furthermore, the /-term in Eq. (22) gives zero contribution within the Glauber theory both for the 
single and double scattering, because for the on-shell pN scattering involved into multistep scattering (28), one has 
(q • k) = 0. For the same reason, the component of the d-term proportional to 0-^(7 vanishes in Eq. (22) too. The 
rest g', g and h terms contribute to the double scattering forward pd-elastic amplitude. 

The double scattering amplitude is given by integration of Eq.(29) over Q = q' [13] 

^ Jj d2<?'M(q, q'; S, zr). (32) 

According to Eq. (11), in order to get the TVPC g amplitude one has to calculate the matrix element of the operator 
given by Eq.(32) at q = 0 over definite initial [p, A > and final \p', A' > spin states: 

5 = J ^ 0lM(q= 0,q';S,cr)|p = = 1 > . (33) 

When considering the double scattering mechanism, in addition to three vectors {k, q, n} defined after Eq. (13) it 
is convenient to introduce two more sets of orthonormal unit vectors {kj,qj,nj} for the first (j = 1) and second 
(/ = 2) collision as it was done in Ref. [13]. At zero scattering angle we have qi = —q 2 = —q', where qi (q 2 ) is the 
transferred momentum in the first (second) collision; rii = — [k x q'], n 2 = [k x q'], ki = k 2 = k -P q'. In the eikonal 
approximation vectors q and q' are orthogonal to k and we assume ki = k 2 = k [13]. The Cartesian projections for 
unit vectors can be written in terms of components of the two-dimensional vector q' = {q^^qy) (OZ axis is directed 
along k) as fii^ = q'y/q', fiiy = -qUd', = 0- 
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FIG. 1: Double scattering mechanism with TVPC (black circle) and T-even P-even (open circle) charge-exchange pn interaction. 

A. g' term 

Non-zero matrix elements of the isospin-operator connected with the g' term in Eq. (22) are 

< n,p\[T X TN]z\p,n >=-i2, < p,n\[T x TN]z\n,p >= i2. (34) 

Therefore, the g' term contributes only to the charge exchange transitions. Two allowed double scattering amplitudes 
with one TVPC and another one T-even P-even pfV-interaction are depicted in Fig. 1. Within the operator formalism 
these two terms can be evaluated in the following way. For pure T-even P-even (TCPC) interactions the transition 
operator for the charge-exchange mechanism of the pd —>■ pd process has a form [29] 


Otvpc — ~ 2 [Mc(q2)Mc(qi)], 


(35) 


where Mc(q) = M„(q) — Mp{q) is the transition operator for the strong charge-exchange pn —>■ np amplitude that 
is equal to the np —>■ pn amplitude. We note that according to Eq. (34), for the TVPC NN-interaction with the g' 
terms the amplitude pn —>■ np differs from the amplitude np —>■ pn by the sign. In order to get the TVPC operator of 
the charge-exchange pd scattering, we make in Eq. (35) the replacement Mc(q) —>■ Mc(q) -I- T'c(q), where the index c 
means either pn —>■ np or np —>■ pn and is the TVPC charge-exchange NN-scattering operator, normalized as 
in Eq. (13) and related to the tp^- operator given by Eq. (22) as 


TpN = 


niN 


4^kpN 


tpN . 


(36) 


Furthermore, we neglect the terms of the second order in Tc as compared to the first order and omit the T-even term 
~ McMc- As a result, the TVPC charge-exchange operator takes the form 




= -i[M„ 


ii<l2)Tpn —^np(Ql) “1“ '^np—¥pn {^.2)Mp 


,(qi)]- 


(37) 


For further evaluation it is convenient to use M^p^pn = Mp^^np = — Mp. Under the sign of the integral over q' 

in Eq. (33) the operator (37) is not changed after the substitution qi o q 2 . In order to find the operators U, Vp, 
V„, and Wij introduced in Eq. (28), it is convenient to add to the right side of Eq. (37) the term O^Yp(j{l O 2) 
and divide the obtained sum by the factor of 2: 


Otvpc — fi + fill 


(38) 


where 


fi = --[(M„(q 2 )Tp„^„p(qi) -b T„p^p„(q 2 )M„(qi)) -b (qi O q 2 )]. 


fll — . [A/p(q 2 )Tpn—>-np(qi) 4” Pnp—>-pn(q 2 )A/p(qi) -b (qi ■f-)- q 2 )]- 


(39) 


Using Eqs. (13), (34) and symmetry in respect of the replacement (qi O q 2 ), we find that Am, Bm, Gn, and Hm 
terms cancel in operators // and fip. 


f - i- 


fii = 


Cn{(T ■ ni)(cr„ - (Tp) ■ m - C(,(cr„ • ni)(crp) • m) -b 
Cp{(T ■ hi){(Tp - cr„) • m - Cp{(Tp ■ ni)(cr„) • m) -b Cpiiini 


(40) 



where 


n = 4v^mArfcpAr. (41) 

Making the sum // + fjj we find the operators U, V, and Wij in Eq. (28) for the g'-term as the following 

[/=^(C; + C'>ini, 

Vp = (Cp - C'„)(<T • ni)ni, V„ = (C„ - C'p)(<T • ni)ni, 

=-lf(^n + Cpnuni,, Wu = + Cpn.n,. (42) 

One can see from Eqs. (42), that V = Vp + V„ = 0. Furthermore, taking into account the relation Vpcrp + V„cr„ = 
VS, we find that for the g'-term the operator Eq. (28) does not depend on the spin of the proton beam cr. As a result 
the transition operator given by Eq. (29) is diagonal in respect of spins of the proton beam. According to Eq. (11), 
it means that the contribution of g' term to the TVPC amplitude g is equal to zero. We emphasize that this result is 
true for the S- and D- components of the deuteron wave function and for all spin terms in the transition amplitude 
(29) allowed in the Glauber formalism. It is easy to find that this result is valid for the nd scattering too. 


B. h and g-terms 

The TVPC interaction corresponding to the h and g terms in Eq. (22) occurs both in the pp- and pn- elastic 
scattering. Following to Ref. [29] (see Eq. (2.7) in it) and [13] we consider the symmetric and antisymmetric 
parts of the operator 

Q+^ = 2 [(^^^(*41) + -^p(qi)(Tpn(q 2 ) + -^n(q 2 )) + (rp„(q 2 ) + Al„(q 2 ))(Tpp(qi) + Mp(qi)], 

Q-^ = 2 [(^pp(qi) + -^p(qi)(^pn(q2) + ^n(ci2)) — (7pn(q2) + -^n(q2))(7pp(qi) + -^p(qi)]- (43) 

According to Refs. [29] and [13] the matrix elements of the antisymmetric operator over the spin S' = 1 states of 

the deuteron are negligible and therefore, we drop this term. In the operator one can neglect terms of the second 
order in the TVPC interaction (~ TppTpn) as compared to the first order Tppf and should omit the pure T-even terms 
~ MpMn- Thus, the double scattering TVPC operator consists of the following four terms 


g(d) ^ hTpp{c[i)Mniq_2) + M„(q 2 )rpp(qi) + Tp„(q2)Mp(qi) + Mp(qi)Tp„(q2)]. 


(44) 


As in the case of g'-term, it is convenient to add to the right side of Eq. (44) the term 0^^(1 o 2) and divide the 
obtained sum by the factor of 2. 

For the h-term we find the operators U, V, and Wij in Eq.(28) for the forward double scattering as 


[/ = 0, Vp = 0, V„ = 0, Wu = 0, 


(45) 


Wij = C'„{q')hpid')[(^^in2iqi] + crqifcijh2i] + Cp{q')hn{q')[(Tkin2jqii + crqikiin2j] 


Wij {Si, Sj} — 


C'nhp + C'„h„ 


{(crki)[(n2S)(qiS) + (qiS)(n2S)] + (crqi)[(kiS)(n2S) + (n2S)(kiS)]} . 


(46) 


In Eq. (28) in this case only the operator Wij depends on the beam proton spin and, therefore, gives non-zero 
contribution to 'g. The matrix elements over the proton spin states are 


{C'^hp + C''/i„)(—-I- iq'p) 


< p' = +^\Wij{Si,Sj}\p = >= 


- S^Syq'^ - q'pS^S, + q'^SyS, 


x^y^z 


n|ni|' 


(47) 
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For the deuteron spin matrix elements of the operator M one has < A' = 0|M|A = 1 >= —{M^x + iMzy )/Finally 
the spin matrix element in the S-wave approximation is the following 


< = 0|M(q = 0,q';S,cr)|/r = -i,A = 1 >= _ ^ +^phn)q 


For the g-term we have got 


[/ = 0, Vp = 0, V„ = 0, Wu = 0, 
Wij{Si,Sj} = ^^" |(n2 • S)([ni x cr] • S) + ([ni x cr] • S)(h2 • S)| 

Furthermore, using the proton spin matrix element 

< = ^l([ni X cr] • S|^ = -i >= {-iq'y + q'x)S^k, 

we find in the S-wave approximation 


< = 0|M(q= 0,Q;S,cr)|/r 


-2,a = i>= 


ik {C'^gp + C'pgn)q' 

72 n 


Si°\q'). 


(48) 


(49) 


(50) 


(51) 


The operator (28) does not depend on the deuteron state. Therefore, the contribution of the D-component does 
not change the factor C'^^hp + C'phn in Eq. (48) and the factor C'^gp + Cpgn in Eq. (51). The D-component can 
contribute due to the last two terms in Eq.(29) with Wij caused by interference with the S-wave and by the w'^- 
term. According to analysis performed in Ref. [13], at energies ~ 100 MeV the D-wave contribution as well as the 
spin-dependent pN- scattering amplitudes are less important than the S-wave contribution and spin-independent pN 
amplitudes. Therefore, we postpone investigation of the D-wave contribution to the next paper. 


C. g amplitude 

Thus, the double scattering mechanism with TVPC interaction from Eq. (22) leads to the following result for the 
g amplitude in the S-wave approximaion 

- / dqq^ S''o\q)\C'^{q){gp - hp) + Cy{q){gr, - hn)]. (52) 

4Wp7r jq 

V. NUMERICAL RESULTS AND DISCUSSION 

The main aim of this study is to analyze the null-test signal a within the Glauber theory. In order to demonstrate 
capability of the Glauber model at energies of the planned GOSY experiment [10] we calculated several spin observables 
of the pd scattering at 135 MeV in comparison with the existing data. The results of our calculations for the unpolarized 
differential cross section, vector Ay and tensor analyzing powers, and spin correlations parameters Cij, Cij^k given 
in Eqs. (31) are in reasonable agreement with the available experimental data and/or Eaddeev calculations [15, 16] 
at 135 MeV and 250 MeV in the forward hemisphere (0cm < 30°). Some of these calculations are shown in Eig. 
2. We also found that Coulomb effects taken into account as explained above improve agreement with the data on 
the non-polarized differential cross section and vector analyzing powers and Ay at these energies at small angles 
9cm < 20° — 30°. The obtained results lead us to the conclusion that the Glauber theory is quite suitable for studying 
the null-test signal for TVPC effects in the pd-scattering because the corresponding signal is not affected by the strong 
background of T-even P-even interactions. 

The previous study of the null-test signal in pd scattering was performed in Ref. [12]. The integrated cross 
section a was calculated within an approach accounting separately for the elastic channel and the deuteron breakup 
pd —>■ pnp estimated within the single scattering approximation. The double scattering mechanism was not considered. 
Extension of the calculation in Ref. [12] to higher energies above the pion threshold is questionable because the mesons 
production is not taken into account in Ref. [12]. On the contrary, our approach based on the optical theorem is 
more general and allows one to overcome these drawbacks. In particular, we find that the a observable is determined 
by the double scattering mechanism. 
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FIG. 2: (Color online) Results of our calculation [14] of the spin observables Cxz,y (a), Cz,x (b), Cy^y (c) and Cx,z (d) for the 
pd elastic scattering in comparison with the data [16] at 135 MeV: without (dashed line) and with the Coulomb interaction 
included (full). 


Our main result obtained within the Glauber theory is formulated by Eq. (52) for the null-test signal. It is worth 
noting that in our approach only the amplitude C'pf appears in Eq.(52). Some other T-even P-even pN amplitudes, 
which were found in Ref. [12] to contribute to the TVPC null-test signal, are absent in Eq. (52). There are two other 
points worth mentioning in relation to Eq. (52). First, the g'-term makes a zero contribution to g and this result is true 
in general case when both the S- and D-components of the deuteron wave function are taken into account. Therefore, 
an exchange by the lightest meson, that is the p-meson, allowed in general case in the TVPC NN interaction [7] and, 
as expected, makes the most important contribution to the TVPC NN interaction, does not contribute to the null-test 
signal a. Contribution of other heavier mesons is usually expected to be less important due to the NN repulsive core 
at short distances between nucleons. A microscopic T-violating optical potential for the nucleon-nucleus interaction 
was derived in Ref. [27] starting from the T-violating p-meson interaction between nucleons. This potential and the 
corresponding coupling constant of the p—meson to the nucleon pp is widely used [ 6 , 7] as a measure of intensity of 
the TVPC effects. However, as we have shown, for the nucleon-deuteron scattering this parameter cannot be applied 
strightforwardly as a scale of the TVPC interactions. 

Strong suppression of the contribution of the p-meson as compared to the axial hi meson was found numerically 
in the Faddeev calculations [25] of the null-test signal for the nd scattering at 100 keV, but no explanation of this 
result was offered. We suppose that the cause for this suppression is the same spin-isospin structure of the scattering 
amplitude which leads to the vanishing p-meson contribution in the Glauber approach. A qualitative explanation of 
the vanishing contribution of the TVPC charge-exchange amplitude for the TCPC terms C'j^ is the following. For 
the strong (TCPC) interaction the charge-exchange amplitude Mpn^npin) appearing in Fig.l,a coincides with the 
charge-exchange amplitude M„p^p„(q) in Fig. l,b. On the contrary, for the TVPC interaction caused by g' term 
the corresponding charge-exchange amplitudes have the opposite signs due to Eqs. (34). The corresponding deuteron 
vertices are the same in Fig. l,a and b. Taking into account the symmetry in respect of the substitution qi -o- q 2 
under the sign of the integral over q' and keeping in mind that spin-dependence of the strong NN scattering amplitudes 
with the C'jq term is identical to that for the TVPC g' term, we find that the double scattering amplitude in Fig.l,a 
differs from that in Fig.l,b only by the sign. Therefore, the sum of these diagrams for the and g' terms is zero 

[33]. 

The second point is connected with the role of the Coulomb interaction in the cross section a. Being a T-even P-even 
interaction the Coulomb pp scattering cannot generate the TVPC amplitude 5 within the single scattering mechanism, 
therefore its contribution to 'g is zero in this approximation. In order to include the Coulomb interaction within the 
double scattering mechanism of the Glauber theory one should replace the pure hadronic T-even P-even pp amplitude 
Mp given in Eq.(13) by the sum Mp + where is the properly normalized Coulomb pp-scattering amplitude 
(16). It is evident from the spin structure of this amplitudes that the Coulomb term is added to the spin-independent 
term Ap, Ap ^ Ap + and double spin terms Bp, Gp, Hp but does not enter into the single spin terms Cp and 
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Cp. However, all amplitudes Ap, Bp, Gp, Hp are excluded from the above derived formula for the TVPC amplitude 
(52) due to the specihc spin structure of Eq. (2). As was noted in Sect. IV B the latter statement is also true if the 
D-wave of the deuteron is taken into account. The only factor in Eq.(52) which contains the Coulomb effects is Cp{q). 
This is because the Coulomb scattering pp amplitude enters into the spin-independent term Ap, which, in turn, enters 
into the amplitude Cp with the relativistic correction factor qj2mp [30] {nip is the nucleon mass) Cp = Cp -I- iAp^^. 
When substituting this amplitude Cp into Eq. (52) and making integration over the transferred momentum q, one 
can see that due to the presence of the factor ^ q^ in the integrand the singularity of the Coulomb amplitude (17) at 
9pp = 0 does not lead to divergence of the g amplitude. 



T [MeV] 

FIG. 3: The calculated energy dependence of the TVPC cross section \a\ for the h-term in units of the constant (j)h = Gh/Gh 
for TUh = 1.17 GeV, A = 2 GeV in /ijv by Eq. (24) with the Coulomb interaction included (dashed line) and excluded (full). 
The dot-dashed curve is obtained at A —>■ oo, -|- ^ m\ without account for the Coulomb effects. 

Energy dependence of a is calculated here for the h term in units of the unknown hNN coupling constant (j)h = 
Gh/Gh with G\/Att = 1.56 and ruh = 1-17 GeV [25]. For the monopole form factor Fhnn we used A = 2 GeV 
[12]. The pN scattering amplitudes Cn and An were calculated using the SAID database [31]. The results obtained 
with and without allowance for Coulomb effects are shown in Fig. 3. One can see that the magnitude of a smoothly 
decreases with increasing beam energy. This is caused be the energy dependence of the strong T-even P-even pN- 
scattering amplitude. The role of the Coulomb interaction in the null-test signal a is rather unimportant in the 
considered interval of energies (Fig.3). At a beam energy 135 MeV we obtain a = 0.039(/)/i mb and (Tq = 78 mb. 
Assuming the accuracy of the experiment allows to measure the ratio ct/cto at the level ~ 10“®, the bound on (f>h can 
be achieved as c/h < 0.002. If we neglect some details of the phenomenological h-exchange and put FhNN = 1 and 
(m| -I- q^)“^ —f ml in Eq. (24), the energy dependence of a corresponds to the dot-dashed curve in Fig. 3. In this 
case we obtain at 135 MeV a = 0.063(j)h mb and, therefore, at the same accuracy of the experiement the bound is 
4>h < 0.0012. The (/-term leads to a very similar energy dependence. Fig. 3 shows that the TVPC cross section a 
increases with decreasing energy. Perhaps, at energies below~ 50 MeV the sensitivity of the experiment [10] to the 
TVPC effects is higher. However, at these low energies the applicability of the Glauber theory to the spin observables 
of the pd scattering is not validated, and therefore rigorous calculations like the Faddeev ones are needed. 

Let us consider possible false-effects in the planned experiment [10]. One source of these effects is connected with 
the non-zero vector polarization of the deuteron Py ^ 0 directed along the vector polarization of the proton beam 
In this case the term aipl/py in Eq. ( 8 ) contributes to the asymmetry corresponding to the cases of PyP^z > 0 
and PyP^z < 0 which is planned to be measured in the TRIG experiment [10]. According to our calculation, at a 
beam energy of 135 MeV the total cross sections are CTq = 78.5 mb, a\ = 3.7 mb, (T 2 = 17.4 mb, and (Xg = —1.1 mb. 
Therefore, the ratio r = (j\/aQ is « 0.05. If the TRIG project is going to measure the ratio Rt = o/a^ with an 
uncertainty about < 10“® (upper limit for Rt), one can End from the obtained ratio r that the vector polarization of 
the deuteron pjj has to be less than « 2 x 10“®. When making this estimation, we assume that the background-to- 
signal ratio is Py crj/a ~ 10“^. The total hadronic polarized cross sections ai (z = 0,1, 2,3) are calculated here using 
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the optical theorem. The Coulomb effects for these observables can be taken into account in the line of Ref. [32] using 
the beam acceptance angle. 


VI. CONCLUSION 

Using the representation of Ref. [18] for the forward elastic pd scattering amplitude and including the phenomeno¬ 
logical TVPC-term in the most general form, we show on the basis of the optical theorem that this term generates 
an extra spin dependent total cross section of the pd scattering. This additional term is zero if the TVPC interaction 
is absent and non-zero only in the presence of the TVPC interaction. Earlier this conclusion was found in differ¬ 
ent representations [II, 12]. Obviously, this null-test signal for T-invariance violation is not affected by the initial 
or/and final state interactions, because it is derived from a genuine pd-scattering amplitude considered beyond the 
perturbation theory. Furthermore, using the Glauber theory we show that (i) the TVPC interaction caused by the 
p-meson exchange does not contribute to the null-test observable a and (ii) the Coulomb interaction does not lead to 
divergence of this observable. Numerical calculation of the energy dependence of a shows that the choice of ~ 100 
MeV made in Ref. [10] is more preferable than ~ 1 GeV. If the cross section a will be measured with the planned 
accuracy [10], the bounds on T-odd coupling constants of the NN-interaction can be achieved as (j)h < ^ 2) x 10“^. 
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